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I. INTRODUCTION 

Mixing processes in fluids play a key role in a wide variety of engineering applications and for 
natural systems such as the ocean and atmosphere. Their theoretical study has been a major focus 
of research, as indicated by the large number of review articles [1-8]. At the smallest scales mixing 
is achieved by molecular diffusion processes, but it may be facilitated greatly by stirring. The 
result of stirring is usually to enhance the effect of molecular diffusion and increase the mixing 
rate [9-13]. Quantitative understanding of the fundamental features of stirring and its influence 
on mixing processes is important for the effective modeling, simulation and design or control of 
these systems. 

The "efficiency" of mixing means different things in different contexts. For example the dis- 
persion of an initial distribution by an imposed flow is a transient problem where the temporal 
approach to the final fully mixed state, rather than the final state itself, is of central interest. 
Consider for definiteness the homogeneous advection-diffusion equation for a passive scalar field 
9(x, t) stirred by a divergence-free velocity field u(x, t), 

r)f) 

+ U -V6 = kA6 (1.1) 

at 

where k is the molecular diffusivity. If this equation is supplied with initial concentration 0(x, 0) 
and applied in an appropriate domain without sources, sinks or scalar flux at the boundaries, then 
the integral of 9 is conserved, so without loss of generality it may be taken to vanish from the start. 
But the L2-norm ||0(-,t)||2, proportional to the scalar variance in finite volume domains, decreases 
with time. Indeed, multiplying (1.1) by 6 and integrating by parts, 

^||0||! = -2s||V0||!, (1.2) 

indicating an inexorable decay of the variance. Efficient mixing in this transient decay problem 
means faster decay of the scalar variance. The mixing efficiency of a particular flow could be defined, 
for example, in terms of its ability to reduce the variance from the initial value to a prescribed 
value within a specific period of time [14]. Because the right-hand side of (1.2) is proportional to 
k, it is evident that molecular processes are ultimately responsible for mixing by this criterion. 
Even though the stirring field does not appear explicitly in (1.2), the conventional intuition is that 
material line stretching in the flow can amplify scalar gradients thereby enhancing the molecular 
mixing rate. Indeed, the velocity's rate-of-strain matrix serves as the local growth rate of the scalar 
gradient field. These issues are of extreme interest for both theory and applications, but in this 
paper we are interested in a distinct scenario where different effects are at work. 
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Mixing a scalar field whose fluctuations are constantly replenished by steady but spatially in- 
homogeneous sources and sinks is a problem with a long history. Early on, Townsend [15, 16] 
was concerned with the effect of turbulence and molecular diffusion on a line-source of temper- 
ature, a heated filament. The spatial localization of the source, imposed by experimental con- 
straints, enhanced the role of molecular diffusivity. Saffman [17] also found that molecular diffusion 
and turbulent diffusion were not simply additive and that higher-order corrections were needed. 
Durbin [18] and Drummond [19] introduced stochastic particle models to turbulence modeling, 
and these allowed more detailed studies of the effect of the source on diffusion. Sawford and 
Hunt [20] pointed out that small sources, such as heated filaments, lead to an explicit dependence 
of the variance on molecular diffusivity. Many refinements to these models followed, see for in- 
stance [21, 22] and the review by Sawford [5]. Chertkov et al. [23-27] and Balkovsky & Fouxon [28] 
treated the case of a random, statistically-steady source. Our goal in the present paper is to 
make the source-dependence of the concentration variance more precise by working directly from 
the advection-diffusion equation, without specifying the underlying turbulent statistics other than 
basic stationarity and homogeneity assumptions. 

When a source of scalar concentration is present, the transient kinetics are of less immediate 
interest and instead the properties of the (statistical) steady state are of greater relevance. As 
will be seen, this sustained steady-state dynamics highlights other features of stirring and mixing 
processes. In comparing the steady-state problem to the transient problem defined by Eq. (1.2), 
it is important to remember that the long-time asymptotic behavior of the decaying problem is 
usually irrelevant to the corresponding long-time behavior of the steady-state problem. This is 
because the continuous replenishing of concentration overwhelms small-amplitude effects observed 
for long-time decay, such as the 'strange eigenmode' [29-40]. 

In this paper we consider the stirring and mixing of a passive scalar sustained by a steady 
source-sink function s(x). Given a prescribed divergence-free velocity field u(x, t) and a molecular 
diffusivity k, the scalar concentration 0(x, t) obeys the inhomogeneous advection-diffusion equation 



supplemented with initial concentration field 6(x, 0). We consider a domain without any net scalar 
flux at the boundaries: the periodic box of size L, i.e., x € T d , the <i-dimensional torus of volume 
L d . The spatial mean of is computed immediately, 



— + U-V6 = kA9 + s(x) 



(1.3) 




(1.4) 
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and deviations from the spatial mean satisfy (1.3) with s(x) replaced by s(x) — L~ d f s d d y. So to 
study the fluctuations we may assume without loss of generality that 0(x,O) and s(x), and thus 
also 9(x,t) have spatial mean zero. 

Fluctuations in the scalar concentration are natually measured in terms of the steady-state 
variance (0 2 ), where we introduce the space-time average. The two averaging operations we use 
are the time average 

- 1 /"* 

F{x) := lim - / F(x,t')dt', (1.5) 

i^oo t J 

assuming as necessary that the limit exists, and the space-time average 

(F) :=^ d J F{x) d d x . (1.6) 

Effective stirring makes the scalar field more spatially uniform, lowering the variance, and this is 
the basic mixing effect that we set out to study. Many investigations have been concerned with 
other statistical properties of the scalar field for this kind of model, such as details of the tails of 
the probability distribution of 9 [2, 6]. While these studies present fascinating mathematical and 
physical issues, in terms of applications they are most likely to be of ultimate use in designing 
closure approximations, i.e., models of the model, in order to accurately estimate bulk measures 
of mixing like variance reduction. In this work we focus directly on the supression of the scalar 
fluctuations as a primary indicator of mixing. 

In terms of the Fourier decomposition of the scalar field, 

&k(t) = ^ f e(x,t)e- ikx d d x (1.7) 



L d . 

where k = (2ir /L)n for n = (m, . . . , n^) € the steady-state variance is 



(o 2 ) - . 

k 

The magnitude k of the wavenumbers naturally index spatial scales, so information about the 
fluctuations at different scales may be obtained by weighting the sums of the Fourier coefficients [41]. 
The simplest indicators of the fluctuations on small and large scales are, respectively, 



<iwi 2 > = ;>> 2 i^ (i.9) 

and 



(|v-^| 2 ) = El4l 2 A 2 (i-io) 
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where the inverse gradient V -1 is defined in Fourier space as multiplication by —ik/k 2 , a well- 
defined operator on these functions with spatial mean zero. For the purposes of this study the 
effectiveness of stirring on mixing at relatively small and large scales will be gauged in terms of 
these norms. Efficient stirring decreases the variances on all scales, although it can be expected 
that any particular stirring may be more effective on some scales than on others. 

It is appropriate here to point out a fundamental and elementary distinction between transient 
and steady-state stirring. For efficient transient mixing the goal is to decrease the scalar variance 
~ ||0(-,t)||2 as quickly as possible by increasing the gradient variance ~ ||V0(-,i)||2 via stirring. 
However in the steady-state problem stirring can only reduce the mean scalar gradient variance — 
and thus the mean rate of variance decay — from its purely diffusive value in the absence of stirring. 

The proof of this (perhaps unexpected) fact is easy. Multiplying (1.3) by and averaging over 
space and time with appropriate integrations by parts produces the well-known balance 

K(\ve\ 2 ) = (6s). (in) 

Inserting the gradient and its inverse on the right-hand side, integrating by parts again, and then 
employing the Cauchy-Schwarz inequality yields 

k(\V9\ 2 } = (0V-V" 1 *) = -<V0-V -1 a) < <|V6'| 2 > 1 / 2 <|V- 1 s|> 1 / 2 . (1.12) 

Note that the long-time solution of (1.3) in the absence of stirring is the steady-state solution 
9q (x) of the diffusion equation with source s(x), 

9 = --A" 1 *, (1.13) 

K 

where the inverse Laplacian is multiplication by — k~ 2 in Fourier space. Solving for the steady-state 
scalar gradient variance in (1.12) and noting that V#o = — k~ 1 V~ 1 s, we conclude that 

<|W| 2 > < <|V0 O | 2 >. (1.14) 

This relationship is uniform in the advecting velocity field u(x, t) implying that no clever stirring 
on any scales can increase the mean scalar gradient variance over its baseline unstirred value. In 
light of this observation we may anticipate that stirring strategies designed to maximize mixing 
efficiency in the sustained source problem are likely to be different from those employed in the 
transient decay scenario. 

The variances (\V P 9 \ 2 ), defined above for p G {—1,0, 1}, are both dimensional and dimension- 
ally distinct quantities. In order to compare them with each other or to compare different physical 
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systems we need sensible nondimensional measures. Hence we define the dimensionless multi-scale 
mixing efficiencies, denoted £ p for p € { — 1, 0, 1}, via 

£ 2 _ (IV^Ql 2 ) a 15) 

S - (| V p(9|2) iLi5j 
where 9q is the steady solution of the unstirred problem defined in (1.13). Effective stirring decreases 
scalar variances relative to those due to diffusion alone, increasing these mixing efficiences. The 
calculation resulting in (1.14) has established that £± > 1. 

Intensifying the stirring often increases the mixing efficiencies and it is important to characterize 
this property in terms of the forcefulness of the flow. The simplest bulk measure of the vigor of 
the velocity field is its mean kinetic energy, or equivalently the rms speed U defined by 

U 2 = (\u\ 2 ). (1.16) 

The nondimensional measure of the strength of the stirring relative to the effect of molecular 
viscosity is the Peclet number Pe that we define as 

Pe = — (1.17) 

K 

using the domain length scale L for simplicity here. It will become apparent that this domain 
length scale may not be the most appropriate one for this purpose; determining the relevant length 
scales is one of the central points of this study. 

For many applications it is useful to know how fluctuations at various length scales may be 
suppressed as functions of k, U, and other features of the problem such as details of the flow 
and source-sink structures. Toward this end it is desirable to know how the multi-scale mixing 
efficiencies depend on the Peclet number, and the notion of "eddy diffusivity" provides a conceptual 
benchmark for this dependence. 

A flow with velocity scale U and "persistence length", "mixing length" or "eddy size" £ that 
characterizes the typical distance a particle travels before changing direction can disperse particles 
diffusively on appropriate space and/or time scales. This suggests that when advection dominates 
molecular diffusion, an effective diffusion with coefficient n c g ~ U£ might replace the advection to 
determine some gross statistical features of the scalar field. 1 If this is so, then the steady-state 
scalar variances are all ~ and the efficiencies are all ~ k c s/k. According to this argument 

U£ I 

^~_ = lxPe. (1.18) 
k L 



1 Indeed, the situation where £ is much smaller than any length scales in the initial data 0(x,0) or the source s(x) 
is the setting for homogenization theory. 
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The linear scaling £ p ~ Pe at high Peclet numbers, which we will refer to as the "classical" 
scaling, provides a baseline reference for the multi-scale mixing efficiencies. Flows that generate 
this classical scaling asymptotically as Pe — > oo produce a truly effective "residual" molecular-like 
diffusion as far as the suppresion of variances at the various spatial scales is concerned, even in 
the singular k — > limit of vanishing molecular diffusivity. And as this discussion suggests, if the 
efficiency scales classically then the prefactor provides a precise prediction of the length scale with 
which an eddy diffusivity might be meaningfully identified. 

The principal purpose of this paper is to determine limits on the mixing efficiencies, i.e., to 
derive a priori bounds on £ p as a function of Pe, and to investigate what sort of flows might realize 
those limits. Upper bounds on £ p (Pe) are of particular interest because they characterize the most 
efficient stirring strategies that can possibly be hoped for. If the bounds are to be useful then they 
should be realizable or approachable, or at the very least indicative of the kind of behavior that is 
possible — such as a scaling like £ p ~ Pe Q at high Peclet numbers. 

Upper bounds on the mixing efficiencies follow from lower bounds on the variances and the first 
study in this direction was apparently by Thiffeault, Doering, and Gibbon [42] who focused on 
estimates and simulations for £q . They adapted an approach that had been used to bound turbulent 
dissipation in the Navier-Stokes equations [43, 44] for application to inhomogeneous advection- 
diffusion equations. For the steady source model of interest here they showed that if (6q) < oo 
then £q< a + 6Pe where the coefficients a and b are homogeneous scale invariant functionals of the 
source s, i.e., invariant under the transformation s(x) — > cs(c'x) for any constants c, and d ^ 0. 

That result showed that very generally the classical scaling is an upper limit to the mixing 
efficiency in this most basic sense. Moreover, the coefficient in the high Peclet number scaling 
£q < bPe puts a limit on any reasonable value for a mixing length: I < bL. It is especially notable 
that this rigorous estimate of t (really of the prefactor b) is uniform in the stirring field u and 
independent of any length scales it exhibits. It also is independent of k. It emerges as a length 
scale in the source-sink distribution, which is seen to play a more important role in the mixing 
process than the conventional eddy diffusion picture anticipates. When the relevant length scale 
in the sustaining source and sink is small then the variance suppression by advection is necessarily 
limited by this no matter what spectrum of scales is present in the stirring process. 

These observations serve as the starting point for this study. Here we carry forward the investi- 
gation of stirring and variance suppression by extending the analysis to multi-scale mixing measures 
while focusing on a broad but specific class of statistically stationary homogeneous and isotropic 
flows. The restriction to this class of flows — a class that includes but is not limited to homoge- 
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neous isotropic high Reynolds number turbulence — allows for the exact solution of some variational 
problems for bounds on the mixing efficiencies. For certain sources and sinks these bounds yield 
anomalous sub-classical exponents for the Pe scaling of some of the £ p . Thus anomalous scaling is 
inevitable for some source-sink distributions. In those cases it cannot be avoided by manipulating 
details of the flow or the spectrum of the stirring. We study the case of a monochromatic source 
s(x) ~ sin(fc s • x) in detail. The estimates are particularly simple in this case and they are sharp: 
we will exhibit a statistically stationary homogeneous and isotropic stirring strategy that saturates 
the upper bounds. 

The rest of this paper is organized as follows. We introduce the relevant class of statistically 
stationary, homogeneous and isotropic flows in Section II, and present some specific examples. 
We formulate variational problems for bounds on the mixing efficiencies in Section III, and derive 
general estimates. In Section IV we evaluate the bounds explicitly for some particular sources and 
sinks and show that classical scaling estimates may be sharp. In Section V we focus on the high-Pe 
behavior of the mixing efficiency bounds for a variety of source-sink distributions, and show that 
anomalous sub-classical scaling is sometimes unavoidable. 

In Section VI we are concerned with the fundamental example of a monochromatic source 
stirred by a flow with a single length scale, the so-called random sine flow, a type of renewing 
flow. Measuring the mixing efficiencies in direct numerical simulations, we find that the large- 
scale efficiency £„i scales (nearly) classically with respect to Pe, like its upper bound, but the 
intermediate and small-scale efficiencies £$ and E\ scale sub-classically, i.e., with powers of Pe less 
than 1. We show that the anomalous exponents for £q and E\ can be deduced from the asymptotic 
analysis of a static flow problem. The concluding Section VII contains a summary of the results 
along with a discussion of open problems and compelling future challenges. Some technical details 
are relegated to appendices. 

II. STATISTICALLY STATIONARY HOMOGENEOUS ISOTROPIC FLOWS 

Our approach to estimating mixing efficiencies is kinematic: the stirring vector field u(x,t) 
is assumed to be given. It could be a solution of the Navier-Stokes equations, or it could be 
a stochastic process with convenient or interesting spectral properties, or it could be a regular 
time-periodic field. The mixing efficiency bounds obtained in this paper will apply so long as a 
few generic statistical conditions are satisfied. Of course not every stirring field satisfying these 
conditions will saturate the bounds, but they are all limited by them. 
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This analysis in this paper is concerned with velocity fields that satisfy the following three 
conditions: 

• The field is divergence free, V • u = 0, everywhere and at all times. 

• The field has finite mean kinetic energy, U 2 = (|w| 2 ) < oo, so that the Peclet number 
Pe = UL/ n is finite. We will also presume more regularity for the velocity whenever necessary 
to carry out formal calculations. This will be apparent in the course of the applications if 
the ultimate estimates depend on other norms of the field. 

• The velocity field is statistically stationary, homogeneous and isotropic. For the purposes 
of this work these qualities are defined by the one- and two-point equal time statistics 
(presuming that these time averages exist pointwise in space) 



«i(av) = (H.l) 

«*(*,■)«,-(*,■) = c i3 { X -y) = (V^) eik<x ~ v) + ^r*« ( IL2 ) 

where C(k) depends only on the magnitude k of the wavenumber k. 

The conditions of finite energy and incompressiblity are familiar and straightforward. In the re- 
mainder of this section we discuss these stationarity, homogeneity and isotropy conditions and their 
relevant implications for the calculation of bounds on the various multi-scale mixing efficiencies. 
We also provide explicit examples, i.e., we describe several flows with these properties that will be 
useful for considerations in subsequent sections. 

First, setting x = y in (II. 2) produces the single point component-by-component covariance 

^ = c y (o) = Ej§ i^- k -vf) + ^r<v < IL3 ) 

Note that because C(k) depends only on the magnitude of the wavenumber, 

fc^O fc^O 

Thus 

Ui(x,-)Uj(x,-) = ^ ££(&)<% = -r-Sij (II.5) 

all k 

where we identify the mean square velocity U 2 = (\u\ 2 ) = J2 C(k). 
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Then if the velocity field is sufficiently regular (or equivalently if C(k) decays sufficiently fast as 
k — > oo) we may also deduce some correlations of the derivatives of u. For example differentiating 
Cij(x — y) in (II. 2) by Xk and setting y = x leads to 



dui(x,-) s-^iC(k) ( kikjk k \ 

because C(k) summed against an odd number of orthogonal components vanishes when the sum 
is absolutely convergent. Differentiating Cij(x — y) by x k and yi, setting y = x, and subsequently 
contracting over k and I gives 



K K fc^O 

where we identify the enstrophy fi 2 := (|Vw| 2 ) = Yl k 2 C(k). In the context of statistical turbulence 
theory the ratio A = U/Cl is (proportional to) the Taylor microscale. 

An example of such a statistically stationary homogeneous isotropic flow is the solution u(x, t) 
of the incompressible Navier-Stokes equations 

du 

— + u-Vu + Vp = uAu + f(x,t) , V-u = 0, (II. 8) 

at 

where v is the kinematic viscosity and f(x,t) is a spatially periodic body force applied to maintain 
a statistical steady state. Of course the forcing would have to be capable of producing the homoge- 
neous and isotropic statistics described above, not an altogether trivial task. This approach would 
be taken to study mixing properties of high-Reynolds number turbulent flows. 

Other flows may be easier and more convenient to implement in simulations or to utilize in the 
analysis of specific models. Some features of the stirring — for instance a stationary energy spectrum 
consistent with developed turbulence — can be realized by specifying the modal amplitudes Wfc(£) 
as mean-zero stochastic processes with k • Uk(t) = and appropriately uncorrelated for distinct 



wavenumbers, to produce any desired energy spectrum C(k) = \uk\ 2 - This is possible when \uk\ 2 
depends only on the magnitude k of the wavenumber k. 

One model of interest is a flow involving just a single wavenumber k ^ 0, the random sine flow, 
a.k.a. the renewing wave flow [29, 30, 42, 45]. (The term 'renewing' or 'renovating' has been used to 
refer to flows that are piecewise-constant in time but change randomly at regular intervals [46-48] . 
Zeldovich [49] introduced a similar single-mode flow, but it was oscillatory rather than renewing 
and thus had poor mixing properties.) In this example the velocity vector field switches periodically 
among steady shearing flows of the form 

u(x) = V2U sin (k-x + tp) (II.9) 
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where k ■ U = 0, \U\ = U and the phase if is selected independently and uniformly from [0, 2ir) 
upon each switch. The switching may be strictly periodic or random in time, and in either case the 
characteristic persistence r is another parameter of this flow. Of course an appropriate selection of 
the directions for k and U must be made in order to realize statistical homogeneity and isotropy 
by the definition used here. 

The kinetically simplest possible example of a statistically homogeneous and isotropic flow is 
one where only (7(0) / 0. At each instant of time the flow is then spatially uniform, a steady wind 
u(x) = U with the direction switching periodically or randomly in time so that the time average of 
u vanishes and the component-component correlation satisfies Uj/u] ~ bij . This wind could sample 
many directions, or as few as 2 x d in the ± directions along an orthogonal set of coordinate axes. 

III. BOUNDS ON THE MIXING EFFICIENCIES 

In this section we derive bounds on the multi-scale mixing efficiencies S v for p = 0, 1, — 1, corre- 
sponding to intermediate, small, and large scales, respectively, for flows that satisfy the statistical 
homogeneity and isotropy conditions in (II. 5) and (II. 6). We first briefly address lower bounds, 
but focus for the most part on upper estimates. The variational formulation and solution for up- 
per bounds on intermediate and large length scales proceeds along similar lines, so we shall treat 
the variance in detail and give a more brisk derivations for the large-scale mixing measure. Two 
different lower estimates on the small-scale variance, corresponding to upper bounds on the mixing 
efficiency at small scales, are derived. One of these small-scale results depends on the spectral 
distribution of energy in the stirring field while the other depends only on the total bulk energy. 

A. Lower bounds on the mixing efficiencies 

Lower bounds on the mixing efficiencies follow from upper bounds on the corresponding vari- 
ances. We already derived a lower bound for the small-scale mixing efficiency £\ in the introduction 
by considering the steady-state variance dissipation-production balance in (1. 11), 

k(\V6\ 2 ) = (s6). (III.l) 

The subsequent result in (1.14), valid for any incompressible flow even without invoking any sta- 
tistical assumptions, is precisely the statement that S\ > 1. 

We expect that this estimate can be sharp in the sense that there may exist a flow field with 
arbitrary Peclet number that may realize it. Certainly this is true in 2d: given the source-sink 
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function s(x) with unstirred scalar distribution 9o(x), just define a flow field with stream function 
ip{x) ~ 6q(x). The streamlines of such a flow are along level sets of 6$ so the flow has no effect. 
Indeed, for this flow u ■ V#o = no matter what the magnitude of U is, so 9q is the stationary 
solution of the advection-diffusion equation for any value of Pe. However, this perfectly "non- 
mixing" flow is not statistically stationary, homogeneous and isotropic, and it is not clear whether 
further constraints derived from the full advection-diffusion equation might be implemented to 
raise this lower estimate for such fluctuating flows. We leave that question for a future study. 

We can follow the same line of reasoning to derive lower estimates on the other mixing efficien- 
cies, although perhaps with less satisfaction. A lower bound on £q requires an upper estimate on 
(9 2 ). Starting from (III.l), recalling that has spatial mean zero and invoking Poincare's inequality 
on the left and Cauchy-Schwarz on the right, 

Afly/SflWI 2 ) 1 ^ < k{\V0\ 2 ) = (s6) < <|V- 1 s| 2 > 1 / 2 <|V^| 2 ) 1 / 2 . (III.2) 

±J 

Hence we deduce the upper estimate on the scalar variance, 

{ ° 2) ^ 4^ (|V ~ 1S|2) ' (IIL3) 



The unstirred variance is 



(el) = ^((A- 1 ,) 2 ), (IH.4) 

rv 



so we have the lower estimate 



,2 = Ml > (( A " lg ) 2 ) = E fc W2vr)- 4 |g fc | 2 

(9 2 ) " L 2 (|V^ S | 2 ) E fe W27r)- 2 |l fe | 2 ' 1 ' 

This lower estimate is strictly positive but because Lk/2ir > 1 it is < 1. Hence it does not rule out 
the existence of flows that might increase scalar variance. Note as well that this estimate depends 
explicitly on the functional "shape" of the source-sink distribution, a feature that we will find for 
many estimates — upper and lower — on the mixing efficiencies. 

The result can in principle be sharpened by solving the variational problem 

(9 2 ) < max{(i? 2 ) | k(|W| 2 ) = (si?)} (III.6) 

where the maximization is performed over all $(x) satisfying the (periodic) boundary conditions 
on the domain. The Euler-Lagrange equation for the maximizer i?*(x) is 

= 2tf* - 2finA$* - fis(x) (III.7) 
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where \x is the Lagrange multiplier enforcing the constraint (III. 1) . In terms of the Fourier coeffi- 
cients the solution of (III. 7) is straightforward, 



but \i is the solution of 



ly ^fc 2 |gfcl 2 =v l^l 2 fIII91 

In general it is difficult to solve (III. 9) for /i, but there is one case where it is easy: if the source is 
"monochromatic", i.e., involves only a single wavenumber of amplitude k s , then fi = —2/nk 2 so 

2 > E^P _ ^Miipf^ _ L (IIU0) 

Hence stirring a monochromatic source can never increase the variance. However if the source-sink 
distribution involves even just two distinct wavenumbers, then the solution of (III. 9) — which must 
be performed numerically — yields a value for /i that produces a lower bound for So that is strictly 
less than 1 [50]. Further details are relegated to Appendix A; this analysis does not prove that 
there actually is some stirring that can increase the scalar variance for such sources, but it leaves 
open the possibility. For the purposes of this study we will settle for this lower bound as far as its 
Peclet number scaling is concerned. That is, 

So > cPe°, (III.ll) 

where the positive coefficient c may depend on the source-sink distribution, and we cannot rule 
out that it might be less than 1. 

A lower bound on the large-scale mixing efficiency follows from the constraint (III. 1) via 
Poincare's and the Cauchy-Schwarz inequalities as well. It follows from (III. 3) that 

^<|V-^) < ™ (,„. 12) 

so that 

1 " <|V^| 2 > - L4 (IV-^I 2 ) E k (Lk/27r)- 2 \s k \ 2 ~ ■ [ ■ 6} 

We note again that for the special case of a monochromatic source, a variational formulation as 
in (III. 6) yields the improved lower bound £_i > 1. It remains an open problem to determine if 
these lower estimates ~ Pe° are sharp for arbitrary sources and sinks stirred by some statistically 
homogeneous and isotropic flow. (We will show that they are sharp in some particular cases.) 
Certainly it will be necessary to use more than just (III.l) to answer this question in general. 
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B. Upper bounds on £ 

This analysis begins by multiplying the advection-diffusion equation (1.3) by a smooth, time- 
independent, spatially periodic "projector function" <p(x) and taking the space-time average and 
integrating by parts to obtain 

= (9{u- V + kA)p) + (ips). (111.14) 

Because this constraint holds for all ip, a lower bound on the variance is 

(6 2 ) > maxmin{(?? 2 ) | = (#(u • Vp + nAip)) + (ips)} (111.15) 

where -&(x,t) varies over all spatially periodic function with unconstrained dependence. 2 This 
min-max variational formulation is equivalent to 



(r) > min{(tf 2 ) | u- W = kA$ + s}. (111.16) 
The multiplier function ip(x) plays the role of a Lagrange multiplier to impose the "Reynolds 



averaged" advection-diffusion equation u ■ Vi? = nA'd + s as a constraint. The formulation of the 
bound as a min-max problem in (III. 15) is, as we will see, convenient for its solution. 

The minimization over # in (III. 15) is equivalent to an application of the Cauchy-Schwarz 
inequality to (III. 14): 

(9 2 ) > max J^'^ N0S = max , (III. 17) 

where we defined the advection-diffusion operator C and its adjoint C\ 

£:=u-V - kA and £ f := -u ■ V - kA . (111.18) 

We explicitly indicate the time average in the denominator of (III. 17) remembering that p is 
time-independent. An important point here is that we can average the time-dependent self-adjoint 
operator LC) before carrying out the maximization over p. 

Maximizing (III. 17) over p is equivalent to minimizing its denominator. Without loss of gen- 
erality, since the functional (III. 17) is homogeneous in <p, we constrain p to have unit projection 
onto the source. Thus we must minimize the functional 



F^(\(p£tf<p-ti((ps-l)), (111.19) 



2 Previously, Thiffeault, Doering & Gibbon [42] derived a bound on this variance without optimizing over ip. 



15 



leading to the Euler-Lagrange equation 



= — = ££t tp - us (111.20) 
where fi is a Lagrange multiplier to enforce the constraint (ps) = 1. The minimizer is then 



J££Q±L 
(«(££t)-',) 



Inserting (III. 21) into (III. 17), we obtain the lower bound 



(0 2 ) > (siCC^s) = (s {k 2 A 2 - V • (mt) + k(2Vm : VV + Au • V)} -1 s>. (111.22) 

Interestingly, this estimate depends only on the mean and equal-point correlation of the stirring. 

Specializing to flows satisfying the assumptions of statistical homogeneity and isotropy in (II. 1) 
and (II. 2) — actually we just use (II. 1) and (II. 5) here — we can carry out the time average in (III. 22), 
yielding 

(0 2 ) > (s(Ctf)-h) = (s {k 2 A 2 - (U*/d)A}-h) = £ K2k J^ 2k2/d - (HI.23) 
Using #o = (— kA)~ 1 s we express this result as an upper bound on the mixing efficiency So- 

£ 2 = m < (s s) _ E fc i^i 2 A 4 fIII241 

(6 2 ) ~ (s(A 2 - (Pe 2 /L 2 d) A)' 1 s) £ fc |s fc | 2 /(fc 4 + fc 2 Pe 2 /L 2 d) 1 ' ; 

where the dimensionless Peclet number Pe = UL/k has been inserted. We observe that like the 
lower estimate on £$ in (III. 5), the upper bound on the mixing efficiency depends on the shape of 
the source-sink distribution. 



C. Upper bounds on E\ 

The gradient variance (|V#| 2 ) can quickly and easily be bounded from below in a similar man- 
ner to the variance. Begin with Eq. (III. 14), integrate by parts, and apply the Cauchy-Schwarz 
inequality to obtain 

(ps) 2 = ((up + KVp) ■ W) 2 < (\uip + kV<p\ 2 )(\V9\ 2 } (111.25) 

so that 

(| W| 2 ) > max- ... . (111.26) 

Nl 1 ' v> (\up + nVp\ 2 ) v ; 
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The right-hand side of (III. 26) is homogeneous in p so we minimize the denominator subject to 
the constraint that (ps) = 1. Under the homogeneity and isotropy assumptions (II. 1) and (II. 5), 
the challenge becomes to evaluate 

min {(K\Vtp\ 2 + U 2 tp 2 } | ((ps) = 1}. (111.27) 

Following a similar development as in Section IIIB, the solution is found to be 

(I V0| 2 ) > (s{-k 2 A + U 2 )-^} (111.28) 

and the mixing efficiency at small scales is bounded according to 

,2 < (s(-A)^s) = E fc l%| 2 A 2 fin 29) 

1 - <.s(-A + Pe 2 AL 2 )-M E k \s k \ 2 /(k 2 + Pe 2 /L*y { " ; 

Like the upper bound for £q in ((III. 24)), this estimate depends on the functional structure of the 
sources and sinks, and on the statistically homogeneous and isotropic stirring only through the 
Peclet number. 

We can improve this upper estimate by avoiding the application of the Cauchy-Schwartz in- 
equality that led to (III. 25). We used that inequality above for expediency, but in fact we expect 
the bound to involve only the gradient (i.e., curl- free) part of the field up + nVp. This can be 
seen by solving the full min-max variational problem 



(|V0| 2 ) > maxmin{(|W| 2 ) | (ps) = ({up + K Vp) ■ Vi?)}. (111.30) 



The minimization is straightforward: 

(|W| 2 ) > max— r (111.31) 

M 1 1 ~ v ((V ■w)(-A- 1 )V ■ w) v ' 

where w = up + nVp. The vector field w can be decomposed as 

w = w- VA^V • w +VA~ 1 V • w (111.32) 



divergence-free curl-free 

where the two components are orthogonal. Only the curl-free part contributes in (III. 31). That is, 
the denominator is just the norm of the curl-free portion: 

((V • w(-A _1 )(V • w)} = ([AA- X (V • u?)](-A -1 )(V • w)) 

= (V(A- X (V • w) ■ V(A- X (V • to)) 
= flVA^V- w\ 2 ) < (\w\ 2 ). (111.33) 
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Hence (III. 31) generally represents an improvement over the expression resulting from application 
of the Cauchy-Schwarz inequality in (III. 26). 

This improved bound depends on the full two-point correlation function of the velocity field: 

((V-w^-A-^V-w) = -L J d d x J d d y V x -w(x,-)G(x-y)Vyw(y,-) 

= j d d x J d d y (-V X V X G) : w(x,-)w(y,-) (111.34) 

where G(x — y) is the Green's function for —A on spatially mean-zero functions with Fourier 
coefficients G(k) = l/(L d k 2 ) for k / 0. Under the assumptions (II. 1) and (II. 2) of statistical 
homogeneity and isotropy, 



Wi(x,-)wj(y,-) = (p(x)f(y)ui(x,-)uj(y,-) + k diip(x)djip(y) 

= Cij(x - y)ip(x)ip(y) + K 2 diip{x)dj(p(y). (111.35) 



In terms of Fourier transformed variables, 

J\2' 



<(v-«o(-A-i)v-«,> = £ ^rrf 1 - w-)i^'i 2 +E + W- ( IIL36 ) 

Thus the mixing efficiency at small scales is actually bounded from above according to 



Si < min Efcl * fc| /k 9 x 



l^fcl 



(111.37) 



This upper bound depends on details of the full spectrum of the stirring velocity field. 

We will not perform the optimization over ip for the general problem here; the implications of 
two-point statistical properties of the stirring on the Peclet number dependence of this bound on £\ 
will be left for future investigations. There is one case, however, where the optimization can easily 
be carried out that shows that the result of (III. 37) may indeed be a quantitative improvement 
over (III. 29). That simple case is when the spectrum of the velocity field is concentrated at k = 0, 
i.e., when the velocity field is at (almost) every moment of time a uniform "wind" in space. For 
such spatially uniform statistically homogeneous and isotropic flows, C(Q) = U 2 while all the other 
C{k) = for k + 0. Then 

f 2 < mi. £fc' gfc l 2 / fc2 V ( u * + kA \a \ 2 - EmMI^ mi ™ 

£l - ™ n E^rf £ + * ) ^ ~ w^w' (IIL38) 

In this case the improvement over (III. 29) is just the extra factor of the spatial dimension d in the 
denominator of the denominator of the denominator in the last term. 
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D. Upper bounds on £_i 

To derive a lower bound on the inverse-gradient variance (IV- 1 ^ 2 ) we begin again with (III. 14), 
insert AA -1 = 1, integrate by parts and apply the Cauchy-Schwarz inequality to obtain 

(<ps) = (V(u ■ Vip + KA(p) ■ VA _1 e) < <|V(«- V¥> + KA^)| 2 )3(|V -1 e| 2 )5. (111.39) 

This gives a lower bound on the inverse-gradient variance 

(IV- 1 ^ 2 ) > max —- ^ - . (111.40) 

Xl 1 ' ~ v (|V(tt- Vp + nAp\ 2 ) v ; 

Recalling that ip is time-independent and restricting attention to statistically homogeneous and 
isotropic flows — assuming as well that the enstrophy VL 2 = (|Vu| 2 ) is finite — the denominator is 

I Vu • V99 + u • VV99 + kVA(/?| 2 = (Q 2 /d)\V^\ 2 + (U 2 /d)(Atp) 2 + k 2 |AVv9| 2 . (111.41) 

Then the optimization over ip in (III. 40) is straightforward. Recalling that 

(IV-^oI 2 ) = (IV-^-^l 2 )/^ 2 (111.42) 

and the definition A = U/Q, we conclude that 

e 2 -i < 



,2 / (IV-^A-^I 2 ) 



(s (-A 3 + (Pe 2 /L 2 d)A 2 - (Pe 2 /A 2 L 2 d)A) ^ s) 

E k \h\ 2 /k* 



(111.43) 



E fc \sk\ 2 /{kS + fe 4 Pe 2 /L 2 d + fe 2 Pe 2 /A 2 L 2 d) ' 
Compare this to (III. 24) for the upper estimate on £q and (III. 29) and (III. 37) for the upper 
estimate on £\. Note that the efficiency depends on the enstropy in the flow, i.e., the flow's 
"shear" or "strain" content, directly through A = U/fl. Interestingly, the (or A) term allows 
for an increase in the bound on the mixing efficiency on large scales via stirring on small scales, 
something that is absent in the bounds on the mixing efficiency at intermediate and small scales. 



IV. SATURATING THE MULTI-SCALE MIXING EFFICIENCY BOUNDS 

The upper bounds on the mixing efficiencies £ p derived in the previous section depend on the 
entire source-sink distribution functions, but on just a few features (U, and for 6-1, A = U/ft) of 
the statistically homogeneous and isotropic stirring field. In this section we show that there is at 
least one combination of sources, sinks and stirring strategies that saturate the upper estimates 
exactly at all Peclet numbers. This establishes that on the highest level of generality the upper 
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bound analysis is absolutely sharp. We also show that there is at least one combination of sources, 
sinks and statistically stationary homogeneous isotropic stirrings that saturate the scaling (if not 
necessarily the prefactor) of the lower estimates on the efficiencies ~ Pe°. 
Consider the simple monochromatic source function 

s(x) = V^Ssmihx!) (IV.l) 

where S is the root mean square amplitude, 2n/k s is the wavelength, and x\ is one of the d 
coordinates. For this monochromatic source the lower bounds on the mixing efficiencies at all 
scales are the same, i.e., for p = —1,0, 1, 

I < £ p (monochromatic sources). (IV.2) 

The upper bounds on £q, £\ and £_i in, respectively, (III. 24), (III. 29) and (III. 38) are generally 
different for monochromatic source-sink distribution functions. But if we restrict attention to 
"uniform wind" flows u(x, t) that are at each instant of time spatially uniform (i.e., X7u(x, t) = 0), 
then f2 = 0, A — ► oo, and we can take advantage of the improved bound for E\ in (III. 38) to see 
that the upper estimates are all the same. For p = —1, 0, 1, 

£ p < \Jl + ~Pe 2 /L 2 k%d (monochromatic sources & uniform winds). (IV. 3) 

In order to show that the upper bounds in (IV. 3) are sharp, we construct a family of statistically 
homogeneous and isotropic flows that approach these limits. The trick is to sustain a uniform wind 
in a given direction for a "long" time so that the steady state is very nearly achieved before 
switching to a uniform wind in another direction. The dynamic transients between the changes in 
flow configurations decay at least at rate nkg , so when the transitions are sufficiently infrequent the 
scalar field is almost always (nearly) in a static configuration. Then we may solve the steady flow 
problem for the scalar field exactly and average the variances over the wind directions to evaluate 
the efficiencies in the limit of slow switches among the directions. This adiabatic averaging method 
can be implemented for other flows, too, as will be considered in Section VI. 

Consider uniform winds of speed U blowing along the "diagonals" , in the 2 d directions given 
by unit vectors ^= (±ei ■ ■ ■ ± e^). Each of these directions is equivalent, so we can solve any 
single problem to evaluate the variances. The steady advection-diffusion equation with source 
s(x) = \/25'sin(/c s xi) and uniform stirring field u(x) = -j= Y^j=i &j i s 
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The solution is of the form 9(x) = Ylj=i F^\ x j) where the functions F^\xj) satisfy the system 
of constant coefficient ODEs 



d 2 F^ U dFW V2S „ 

H sm(k s xi/L) = 0, 



dx\ y/dn dxi k 

cPF® U dF® 



(IV.5) 

0, for 2 < j < d, 



dx 2 yfdn dx n 

with periodic boundary conditions F^(0) = F^{L). Recalling that 9 has, without loss of gener 
ality, spatial mean zero, the solution is 



nk 2 + U 2 /dK 
= 0, for 2<j<d. 

The variance is thus 



sin(fc s xi) cos(/c s xi) 

Kk s V d 



(IV.6) 



( * 2) = ^iTwm - <IV - 7) 

and because the scalar field is monochromatic, the small-scale and large-scale variances simply 
satisfy k 2 (\V~ 1 9\ 2 } = (\V9\ 2 )/k 2 = (9 2 ). The multi-scale mixing efficiencies are all then 




precisely as in (IV. 3). 

On one hand this result is fairly intuitive: the most efficient way to reduce the variance (on any 
length scale) is to direct the flow from source regions directly toward the closest convenient sink 
regions and the likewise from sinks toward sources — i/this can be accomplished effectively given 
the constraints of incompressibility and statistical homogeneity and isotropy. This can be done 
simply for these monochromatic flows on the torus, and we have discovered that such flows actually 
suppress the scalar variance at all scales as well as is possible for any statistically homogeneous 
and isotropic flow field. (Note: this example was inspired by Plasting & Young [51] who showed 
that the steady direct flow across the source saturates the variance bound for single-wavenumber 
sources among all flows regardless of any statistical considerations.) On the other hand this type of 
sweeping flow is somewhat pathological in the sense that it simply transports the source onto the 
sink without "stirring" and "mixing" by the usual meanings of those words. In some cases, such as 
transient mixing as discussed in the introduction, the role of effective stirring is to stretch material 
lines and amplify gradients of the passive scalar to accelerate the action of molecular diffusion to 
dissipate variance on small scales. 
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It is not presently clear whether the upper bounds in (III. 24), (III. 37) and (III. 38) can be sat- 
urated for more general sources and sinks. Moreover, this particular saturation result depends 
explicitly on the geometry of the domain. In geophysical and astrophysical applications, for exam- 
ple, it is natural to consider the domain to be the surface of a sphere. While such a domain admits 
similar concepts of statistical homogeneity and isotropy for the flow and "monochromaticity" for 
the source-sink distribution (eigenfunctions of the Laplacian on the sphere), there is no analogous 
uniform sweeping flow as there is on the torus. This makes it clear that the optimal stirrer is a 
function of both the source shape and the domain. Formulating the optimization problem for the 
best stirring field for a given source-sink distribution remains a problem for future investigation. 

We close this section by pointing out that a similar statistically homogeneous and isotropic 
uniform wind can also be arranged to saturate the high-Pe scaling of the lower bounds on the 
mixing efficiencies. For the monchromatic source (IV. 1), consider uniform winds oriented along 
the coordinate axes in the 2d directions ±ej, switching only after blowing each way for a long 
enough time for the transients to be negligible for the time averages. The wind reduces the scalar 
variances below the variance in 6q only when it blows in the ±ei directions, which occurs 1/d of 
the time. Averaging the scalar variance adiabatically over the wind directions yields the mixing 
efficiencies 



/ 1 + Pe 2 /fc 2 L 2 
P ~V l + (d-l)Pe 2 /dklV { ' 



These efficiencies are monotonically increasing in Pe but bounded according to 



1 < £ P < lim J 1 + Pe2/fc | L2 - = A OIxPe . (IV.IO) 

" p p™ Y 1 + (d - l)Pe 2 /dk 2 s L 2 V d - 1 V ; 

V. HIGH-Pe BEHAVIOR OF THE MIXING EFFICIENCY BOUNDS 

Here we examine the high Peclet number behavior of the upper bounds on the multi-scale mixing 
efficiencies for statistically homogeneous and isotropic flows. As a point of reference we recall that if 
there is an effective eddy diffusion associated with the flow field, so that the stirring suppresses the 
scalar variances in the manner of enhanced molecular diffusion in the form of an equivalent (eddy) 
diffusion ~ U£, then the efficiencies would scale "classically", as £ p ~ Pe 1 , as Pe — > oo. Classical 
scaling allows for the precise identification of equivalent diffusivities k^*^ := k£ p =: U l p that serve 
as the definition of associated "mixing lengths" l p that are independent of the magnitude of U 
and k. The upper limits on £q in (III. 24), £\ in (III. 29) and £-\ (III. 43), reproduced immediately 
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below for reference, all depend on the full structure of the source-sink distribution: 



/ |sfcl 2 A 2 



c < / Sfc l^fcl 2 /^ 4 , v 9 , 



E fc |*fc| 2 /(* 6 + A: 4 Pe 2 /L 2 (i + k 2 Pe 2 /\ 2 L 2 d) ' 



The high-Pe behavior, though, can be discerned from just a few features of the sources and sinks. 

A. Square-integrable sources and sinks 

Consider the case where the source-sink distribution function is square integrable, s(x) G 
L 2 (T d ), so that the Fourier coefficients are square summable: 

^2\s k \ 2 < oo. (V.4) 
k 

Then the Pe — > oo asymptotic behaviors in (V.l), (V.2) and (V.3) are elementary to evaluate: 



/V lsJ 2 /P / max ) 



/ V I £, |2 /max) 



There are three significant features of these scaling bounds worth noting. 

The first point is that these upper estimates all scale classically, allowing us to identify the 
largest possible values for meaningful mixing lengths that we have labeled ^ max ' ) . 

The second remarkable fact is that the largest possible mixing lengths relevant to small and 
intermediate scale fluctuations do not depend on the flow field, but rather only on the source- 
sink distribution. Indeed, £^ max - ) < £^ ) max - ) < L/2ir are particular length scales in the source-sink 
function that have nothing to do with the stirring field or any length scales in the flow. This is in 
direct conflict with the notion that a mixing length should be a characteristic persistence length or 
eddy size in the velocity vector field. Rather, these mixing efficiencies are ultimately limited by the 
structure of the sources and sinks. When ^ max - ) and ^[ | max ' ) are "small" (i.e., when 4 max) < ^ max) « 
L) then the mixing efficiencies are also "small" and not subject to any further improvement by any 
sort of clever stirring designed to enhance the variance reduction at intermediate and small scales. 
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The bound on the mixing length £_i for large-scale variance reduction, however, does depend 
on the spectrum of length scales in the flow through (the Taylor microscale) A. It is interesting 
to note that the bound on £^ ax ^ is an increasing function of A -1 , allowing for the possibility that 
small-scale stirring could enhance large scale mixing. That is, this analysis does not preclude 
small-scale stirring from suppressing large-scale fluctuations in ways that it cannot decrease the 
variance at intermediate and small scales. 

We note as well that the improved upper bound on the small-scale mixing efficiency £\ in 
(III. 37) generally produces an even smaller estimate for ^ max - ) that depends on the spectrum, i.e., 
the magnitude and distribution, of the length scales in the flow. In cases where the shortest length 
scale (call it £ s ) in the source is much longer than the longest length scale (call it £ u ) in the stirring 
field, a thoughtful examination of (III. 37) suggests that ^ max ) ~ £ u <^ £ s . 

The third point worth noting is that the upper estimates for all three of the £ p scale the same, 
~ Pe 1 as Pe — > oo. So far all the examples we have considered share this property, but in the next 
subsection we will see that this is not a general feature of the multi-scale mixing efficiency bounds. 

B. Measure-valued source-sink distributions 

If s(x) ^ L 2 (T d ) then the sum in the denominator of (V.5) diverges and the ratio defining 
the prefactor (^ max ^/L) to the Pe 1 scaling could vanish. This would violate the lower bound 
£i > 1, so the high-Pe asymptotic analysis of the ratios of sums must be revisited. This issue 
arises for measure- valued sources and sinks, i.e., when the distribution involves singular objects 
like 5-functions. Then "anomalous" sub-classical Peclet number scalings for some efficiencies are 
inevitable. 

Consider the extreme cases where the Fourier coefficients of s(x) obey 

\s(k)\ = 0(1) as | As | -► oo. (V.8) 

In this case, in spatial dimensions d = 2 and 3, sums in both (V.5) and (V.6) diverge so the high- 
Pe behavior of the bounds on £\ and £q must be re-evaluted directly from (V.l) and (V.2). The 
classical high-Pe scaling (V.7) for the bound on £_i from (V.3) remains the same in d = 2 and 3. 

To discern the high Peclet number behavior of the expressions in (V.l) and (V.2) when |s(fc)| ~ 
C = 0(1) as |fc| — ► oo, we analyze integral approximations to the sums. This is justified because 
the major contribution to these diverging sums comes from the high-fc end where the discreteness 
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of the wavenumbers is negligible. For example, the denominator for the bound on £q in (V.2) is 



E 



\ s k\ 



k 4 + k 2 Pe 2 /L 2 d 



_L\ d f°° C k d ~ 3 dk 
>J d J 27T/L k 2 + Pe 2 /L 2 d 



(V.9) 



where S2 = 2n and S3 = Ait. 
In d = 2 this is 



E 



l s fc| 



^ fe 4 + k 2 Pe 2 /L 2 d 
so the upper bound on £0 in d = 2 is 

Pe 



So < 



8tt 2 log 



1 + 



Pe 2 



1? I? 

— C —* log 

2vr Pe 2 



Pe 

(logPe) 1 / 2 



1 + 



Pe 2 

8vr 2 



(V.10) 



as Pe — > 00. 



(V.ll) 



This upper bound exhibits a logarithmic correction to classical scaling. The significance of this 
is perhaps best appreciated as the absence of any residual variance suppression in the limit of 
vanishing molecular diffusivity. That is, the largest possible effective diffusivity defined by the 
bulk variance suppression k c q < k£ ~ OdlogKl-V 2 ) as k — > 0. It is worthwhile stressing 
that this result does not depend on — and cannot be circumvented by manipulating — any further 
details of the statistically homogeneous and isotropic stirring field. 

The deviation from classical scaling is more dramatic in d = 3 where 

12 

.2- J 



E 



\ s k\ 



^ fc 4 + k 2 Pe 2 /L 2 d 
Then the upper bound on £q is 



L V An ^ L 



7T 



arctan 



I Pe I 



(V.12) 



So < 



Pe 



f " arctan (^) 



Pe 1/2 as Pe ^ 00. 



(V.13) 



This upper bound exhibits strictly sub-classical scaling at high Peclet numbers. 

Both the numerator and the denominator diverge in the bound in (V.l) for £\ for such measure- 
valued source-sink distributions, so those bounds must be evaluated as the limit of ratios for a 
sequence of appropriately mollified sources and sinks. This is straightforward if we simply truncate 
the Fourier series for s(x) at small length scale l s <tC.L and study the ideal case where \s(k)\ = C 
for 2tt/L < \k\ < 2tt/£ s and \s(k)\ = for \k\ > 2ir/£ s , and then take the limit i s /L -> 0. 

Again we analyze integral approximations to the sums; the numerator in (V.l) is 

2vr/ J 27 r/L 



E 



|£fc| 

k 2 



k 2 



(V.14) 
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and the denominator is 



E 



l«fc| 2 



k 2 + PeVL 2 



C k d ~ l dk 
2n/L k 2 + Pe 2 /L 2 



(V.15) 



In d = 2 spatial dimensions this implies the bound 



Si < \ 1°S 



L 2 



L*2 



log 



Pe 2 + 4tt 2 



L 2 



as 



0. 



(V.16) 



Pe 2 + 4tt 2 

Hence there can be no reduction of the small-scale variance beyond that due to molecular diffusion 
for such singular source-sink distributions, no matter what flow strategy is adopted or how energetic 
the stirring is. 

It is interesting to note as well that if £ s is small but finite and Pe S> i- S> 1, then the bound 
scales classically again. That is, as Pe — > oo, 



£l < A/log 



L 2 



log 



Pe 2 + 4vr 2 ^ 
Pe 2 + 4vr 2 



2vr v / i 2 - i 2 s 



Pe 



4, 



2vrL 



Pe. 



(V.17) 



This can be interpreted as the statement that the largest possible value for the mixing length 

^(max) _ Q^g^ j n sucn case s. 

For such distributions in d = 3 spatial dimensions, 

-V 2 



£ i < i 



Pc 



- 1 



arctan 



2vrL 
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arctan 



2tt 
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as 



0. 



(V.18) 



Again, there can be no reduction of the small-scale variance beyond that due to molecular diffusion 
for this kind of singular source-sink distributions, no matter what stirring strategy is adopted. 
Furthermore, when £ s is finite and Pe S> j- 3> 1, the bound scales classically again: as Pe — > oo, 

v -1/2 



£i < i - 



Pe 



arctan 



/2vrL\ 



/2tt\ 
arctan — — 



3 f-1 



4- 2 £ - 1 



Pe 



vg4 

2vrL 



Pe. 



(V.19) 

Not unexpectedly, the largest possible value for a mixing length for this kind of mollified distribution 

is 4 max ) = o{l a ). 

Classical scaling for £q is also eventually recovered for these regularized sources and sinks when 
Pe 3> L/£ s , albeit with very different relationships between the smallest source-sink length scale 
£ q and the maximal mixing length 4, max) . When \s(k)\ = C > for 2vr/L < \k\ < 2ir/£ s and 
\s(k)\ = for \k\ > 2tt/£ s the numerator of the bound in (V.2) is 

(L_\ d s I 1 C k''-hlk 
2vr/ J 27 r/L 



E 



Igfel 
fc 4 



fc 4 



(V.20) 



26 



and the denominator is 



|sfc| 2 



^ A; 4 + k 2 Pe 2 /L 2 d 
For d = 2, as Pe — > oo, 



L 
2^ 
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So in spatial dimension ei = 2 the largest possible mixing length ~ L/ y/\og[L/l s 

length scale vanishes as i s — > 0, but much slower than the actual smallest scale 0(l s ). 
In cZ = 3, 



(V.22) 



This 



Pe 
<W2 ' 



■p as Pe 



oo 



(V.23) 



so £ max ^ = C(^y 2 ). This vanishes faster than in 2-d as £ s — > 0, but still much slower than 

A log-log plot of the bound on £q vs Pe for these cutoff source-sink distributions is shown 
in Figure 1 for the 3-d case with £ S /L = 1CT 8 . The anomalous Pe 1//2 scaling persists over the 
range 1 < Pe < L/£ s and classical scaling ~ Pe 1 takes over at higher Peclet numbers with a 
"small" prefactor ~ ^£ S /L. This small prefactor is a quantitative indication of the difficulty of 
any statistically homogeneous and isotropic flow field to efficiently suppress the scalar variance — 
beyond the suppression achieved by molecular diffusion alone — in the presence of small-scale scalar 
sources and sinks. 

The scalings for the mixing efficiencies with a ^-function source are apparently not all realized 
by a simple uniform wind as was the case for smooth monochromatic sources. We have performed 
the adiabatic approximation for the bulk variance in the case of a uniform wind blowing past a 
5- function source to estimate the mixing efficiencies (see Appendix B). While £\ is necessarily equal 
to 1 in both 2-d and 3-d, we find that the uniform wind gives So ~ a/Pc- (rather than the bound 
Pe/Vlog Pe) in 2-d, and So ~ VPe/^\og Pe (significantly closer to the bound \/Pe) in 3-d. For 
the large-scale efficiency, however, the uniform wind produces S-\ ~ y/Pe in both 2-d and 3-<i, far 
below its classically scaling upper bound. 

A similar integral analysis can be carried out for "fractal" source-sink distributions where \§k\ ~ 
|fc|~ 7 with 7 > 0. The scalings of the mixing efficiency bounds are summarized in Table I. For 
7 > d/2 all the bounds scale classically since then s(x) G L 2 (T d ). The case 7 = is the <5- 
function studied above. We observe that the bounds for S\ and So can exhibit anomalous scaling 
with exponents depending on the fractal nature of the source-sink distribution as characterized by 
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Pc 

FIG. 1: So as a function of Pe for cutoff source-sink distributions, after [52]. 

TABLE I: Scalings of the bound on the mixing efficiency £ p as functions of the source roughness exponent 7 
of the source in two and three dimensions. 



d = 2 p=l p = p = -l 

7 = 1 Pe/^ogPe) 1 / 2 Pe 

< 7 < 1 Pe 7 Pe Pe 

7=1 Pe/^ogPc) 1 / 2 Pe Pe 

7 > 1 Pe Pe Pc 

d = 3 

7 = 1 Pe 1/2 Pc 

0<7<l/2 1 Pe 7+1/2 Pe 

7=1/2 1 Pe/^ogPe) 1 / 2 Pe 

l/2< 7 <3/2 Pe 7 ~ 1/2 Pe Pc 

7 = 3/2 Pe/tlogPe) 1 / 2 Pe Pc 

7 > 3/2 Pe Pe Pe 



7 [52]. Of course if the fractal scaling of the source-sink distribution persists over a broad but finite 
range of wavenumbers, say for 2tt/L < k < 2ir/£ s , then classical scaling for the efficiency estimates 
will again appear for Pe » L/£ 8 . 
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VI. A SINGLE-SCALE SOURCE STIRRED BY A SINGLE-SCALE FLOW 

In the previous sections we have seen that the classical high-Peclet number scalings for the multi- 
scale mixing efficiencies, £ p ~ Pe, are generally upper bounds that may in fact be saturated for 
particular source-sink distributions stirred by certain statistically homogeneous and isotropic flows. 
It was shown in Section V that for some distribution- valued sources and sinks, some of the mixing 
efficiencies necessarily scale anomalously, i.e., £ p < Pe ap with some a p < 1. For those examples the 
bounds on the mixing efficiencies at the different scales scale differently: < a\ < ao < a_i = 1. 
It is not presently known if the anomalously scaling bounds for distribution-valued sources and 
sinks are sharp or, if they are, what kinds of statistically homogeneous and isotropic flows might 
be required to realize them. Those questions remain open, but in this section we settle the issue of 
the possibility of realizing distinct scaling exponents for the mixing efficiencies on different length 
scales for a smooth source-sink distribution. 

The random sine flow, a.k.a. the renewing wave flow characterized by a single length scale, is a 
popular and convenient test bed for studying — both analytically and via simulation — a wide variety 
of stirring and mixing phenomena [29, 30, 42, 45]. Here we consider the simplest single-scale source 

s(x) = >/2S'siii[27r(x + y)/L] (VI.l) 

stirred by a random sine flow that switches between 

u(x) = iV2Usm[2iry/L + 0] (VI.2) 

and 

u(x) = jV2Usm[27tx/L + <j>] (VI.3) 

at time intervals of length r/2, where the phase (f) is chosen randomly and uniformly from [0, 2ir) 
at each switch. This is a statistically homogeneous and isotropic flow, but it is not the maximally 
efficient flow for this source-sink distribution; the optimal flow is the spatially uniform flow con- 
structed in Section IV. We have performed high resolution direct numerical simulations of this 
system in d = 2 spatial dimensions that reveal that the multi-scale mixing efficiencies may have 
distinct scaling exponents. Additionally, for this example we are able to compute the exponents 
theoretically in the context of the adiabatic averaging approximation utilized in Section IV. 

Figure 2 shows the results of the direct numerical simulation for p = 1, 0, —1, each plotted along 
with the classically-scaling upper bounds from Section III. The Peclet number in these simulations 




FIG. 2: Mixing efficiencies £ p as a function of Pe for (a) p = 1, (b) p = 0, and (c) p = —1 for the random 
sine flow with source ~ sin k s (x + y) . The solid lines are the upper bounds for this source from Section III 
and the dashed lines are the data from the direct numerical simulations with U, L and r fixed. 



was varied by holding U = I fixed in a box of side length L = 1 and changing the molecular 
diffusivity k. The computations were carried out with switching time r = 1. The details of the 
numerical simulations can be found in [42]. Writing £ p ~ Pe Qp as Pe — > oo, it is clear that the data 
are consistent with ao < a\ < 1 and a_i f» 1. 

These simulations establish two important facts: (1) that the mixing efficiencies at different 
length scales generally scale differently at high Peclet numbers, and (2) that anomalous subclassical 
scaling can easily be realized by simple and "reasonable" flows. 

We can evaluate the observed scaling exponents a p theoretically for this particular system by 
appealing to a quasi-static adiabatic approximation introduced in Section IV for the optimally 
mixing uniform flow. That is, we solve the time-independent problem for the scalar distribution 
under the influence of steady flows of the form (VI. 2) and (VI. 3), compute the variances, and 
average over the flow configurations. The physical justification for this approximation comes from 
examining the results of the direct numerical simulations: it is observed that the relevant steady- 
state configuration is quickly approached in the time intervals between the switches of the stirring 
field. If we make the switches at longer and longer intervals, i.e., if we increase r, the scalar field is 
distributed (nearly) according to the static configurations for most of the time, and the steady-state 
variances dominates the time averages. 3 

The simplest problem in this category is when the velocity field is oriented parallel to the 
gradient of the source, so we consider the steady advection-diffusion equation 

V2U sin k u y d x 9 = n{d 2 x + d 2 )6 + \/2Ssmk s x (VI.4) 



3 This approach is similar in spirit to "rapid distortion theory" in turbulence [53, 54]. 
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where we allow for different length scales in the stirring and the source with the non-dimensional 
number r = k u /k s gauging the relative amount of shear in the flow. This is not exactly the static 
problem corresponding to the dynamic simulations where the flow is always at a 45° angle from 
the source-sink alignment, but it turns out that the scaling exponents a p are the same for that 
case and for the 90° alignment in (VI. 4). For illustrative purposes it is convenient to explain the 
most elementary case (VI. 4) in detail here and relegate details of the 45° problem to Appendix C. 
The solution to (VI. 4) takes the form 

0(x) = f(y) s'm(k s x) + g(y) cos(/c s x) (VI.5) 

where the functions / and g are periodic on y G [— ir/k u , ir/k u \ and satisfy the system of ODEs 

-V2Uk s sm(k u y)g(y) = K -k 2 s + f{y) + y/2S, (VI.6a) 

V2Uk s S m(k u y)f(y) = « \-k 2 s + -^1 g(y) . (VI.6b) 



From (VI. 6a) we deduce that g(y) is an odd function of y and f(y) is an even function of y, 
compatible with (VI. 6b) which says that / and g have opposite parity. We can thus infer boundary 
conditions on the reduced domain [0, 7r/k u ]: 



g(0) = = g(w/k u ), f'(0) = = f\n/k u ). 



(VI.7) 



We are interested in the high-Pe behavior of the solution to (VI. 4). In Figure 3 we show a 
grayscale plot of the scalar field at several values of the Peclet number. As is evident, internal 
layers develop along the lines of maximum shear around y = nir/k u for n = 0, ±1, ±2, .... Away 
from these lines the flow directly blows source regions onto sink regions and vice versa which, as 
we have seen in Section IV, is the most efficient stirring to reduce variance at all scales. Even 
though the regions of high shear effectively stretch material lines, the mixing process is apparently 
frustrated by the constant replishment of the scalar variation by the steady sources and sinks and 
the variance is dominated by fluctuations concentrated in shear layers. 

We develop an asymptotic singular perturbation internal layer approach for the limit Pe 3> 1 
with r = 0(1) fixed. Upon rescaling y = k s y, f = fUk s /S, g = gUk s /S, r = k u /k s , and 
introducing the slightly re-scaled Peclet number Ve = V2U /nk s = \^2Pe/k s L in Eqs. (VI. 6), we 
obtain the non-dimensional ODEs 



1 

Pe 



1 + — 

dy 2 



f(y) + 1 = -sm(ry)g(y), 



1 

Pe 



-1 + 



dy 2 



g{y) = sm(ry)f(y). 



(VI.8a) 
(VI.8b) 
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(a) 



0.2 0.4 0.6 0.8 



(b) 



0.2 0.4 0.6 0.8 



(c) 



0.2 0.4 0.6 0.8 



FIG. 3: The scalar field stirred by a steady sine flow for (a) Pe = 0, (b) Pe = 100, (c) Pe = 1000. 



Proceeding as usual we construct inner and outer solutions. The outer solution, valid away from 
the internal layers, is obtained by expanding in powers of "Pe -1 : 



/outer ^ ^ fni 



n=0 



5outcr — ^ 9n 
n=0 



To leading order the solution to (VI. 8) in the outer region is 



outer 



> pouter 



sin(ry) 

For the internal layer we expand in a small parameter e as 



E en /« 



E ^ 



n=-l 



dinner — / c :in ■ 
n=-l 



(VI.9) 



(VI.10) 



(VLll) 



so that both /i nne r and gi nne r are C(e _1 ) as e — » 0. The internal layer scaling is determined by a 
dominant balance argument: we choose e = "Pe -1 / 3 and rescale y = erj to achieve a self-consistent 
scaling of the leading order terms. When these scalings and (VI. 11) are inserted into (VI. 8), the 
problem to solve at leading order 0(l/e) is 



d 2 f-i 



+ rrjg-i + 1 = 0, 



rrff-x = 0. 



(VI.12) 



drf drj 2 
Then letting £ = r 1 / 3 ^, F = r 2 / 3 f_i, and G = r 2//3 <7_i, (VI.12) simplifies to the inner layer 
equations 



F +£G+1 = 0, G -£F = 0, 



(VI.13) 



with boundary conditions F'(0) = and G(0) = 0. The other boundary conditions come from the 
requirement of matching to the outer solution (VI. 10): F(£) — ► and G(£) — ► — 1/£ as £ — > oo. 
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FIG. 4: Comparison of the direct numerical solution (solid) and the internal layer solution (dashed) for 
Pc = 1000. 

The system (VI. 13) can be cast as a complex Airy equation for F + iG, but we resort instead to a 
numerical shooting method to obtain the solution. 

The solution of the internal layer equations is obtained by shooting backward (which is the more 
stable evolution direction) from £ — > 00. The large-£ asymptotic behavior deduced from (VI. 13) is 



A + G 

-£4 + £l0> 



1 b 



(VI. 14) 



where in practice a and b are tuned numerically to produce a solution that satisfies the boundary 
conditions at £ = 0. This is easily accomplished, and in Figure 4 we compare the internal layer 
solution against the exact (numerical) solution of the full advection-diffusion equation (VI. 4) at 
Pe = 1000 where the small parameter e ~ .2. This agreement confirms that the leading terms 



/inner — 



-2/3 



no, 



-2/3 



Sir: 



G(0- 



(VI.15) 



do indeed accurately capture the asymptotic behavior. 

The uniform asymptotic solution to the coupled ODEs are, to leading order, the composites of 
the inner and outer solutions. Recovering all the scalings and letting 5 = e/r 1 / 3 /c s , we define 



5c ° mp(y) = ^-^X5inne r (y)X5oute r (y) = ^ ^ Jjg-y G (f) 



(VI. 16a) 
(VI. 16b) 



Uk s k u 5 sm(k u y) 

Armed with the approximate solutions (VI. 16) we can compute the multi-scale mixing measures 
(|V p 6f ) for p = 0,1,-1. 
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The variance is 

i ,,,, / r /Aku ., . r /4ku 



2 ) = ^ ((/ 2 ) + (5 2 )) = = (jf /( y ) 2 ^ + / o </G/) 2 d*/J (VI.17) 



where we use the symmetry of the solution to carry out the integral over only a quarter-period. 
Letting 77 = y/S and replacing / and g by f comp and ff comp , 



^-likUL F "^ + l c( " )2 TO) d " 1 ■ (VU8) 



Here we are interested in the scaling as <5 — > so we are justified in replacing the upper limit 
of the integral of F(rj) 2 by infinity Slightly more care must be taken with the integral involving 
G(i]) 2 : because k 2 r] 2 S 2 / sm 2 (k u r]5) is uniformly bounded by and it/2 we can appeal to Lebesgue's 
dominated convergence theorem to deduce 

1 <?2 1 / roc f-oo \ 

^ ~ f <"» 2d " + l °<"» 2d ")' (VU9) 

The unstirred variance is (0q) = S 2 /n 2 kg, so the mixing efficiency is 

/~ / fOQ fOO \ — 1/2 

- Mf F{r ] ) 2 d 11 + G( V ) 2 d V ) r^Ve^. (VI.20) 

Figure 5 shows £0 as a function of Pe from the numerical solution of the steady advection- 
diffusion equation (VI. 4). The data accurately match the Pe 5 / 6 scaling and using the prefactor 
calculated from (VI.20) there is only a 5% discrepancy between the internal layer solution and that 
calculated from the direct numerical simulation at Pe = 1000 with r = 1. Remarkably, this Pe 5//6 
scaling also fits the data for the time-dependent stirring of the tilted source shown in Figure 2. 

The gradient variance in the steady solution is given by 

<|W| 2 > = f [(f) + (9 2 )] + \ [<(/') 2 > + <G/) 2 >] • (VI.21) 

The (f 2 ) and (g 2 ) terms were computed above in the context of £q, so we focus on the derivative 
terms. Inserting the composite asymptotic solution (VI. 16) we obtain 

\((f>)2\ + ((„i) 2 \] „l (— — Lp'r„AfiV 

/ N 2\ 




1 7 G'( y /6)^f- + JL-L G {v/S) 



k u y 



Uk s k u 5 2 sin(py) Uk s k u 5 [sm(k u y)_ 

The last term is smaller (0(<5 -1 )) than the immediately preceeding term (0(5~ 2 )) and may thus 
be neglected. The first two terms in (VI.21) are also 0(5^ 1 ), so they may be neglected as well. 
Hence 
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FIG. 5: Mixing efficiencies £ p for p = 1 (x), p = (*) and p = — 1 (+) for the steady sine flow. The solid 
lines are the scaling predictions of the asymptotic analysis and the discrete data are from the numerical 
solution of the steady advection-diffusion equation (VI. 4). 

and by the same arguments as for the variance above we deduce 

Using (|V#o| 2 ) = S 2 /K 2 kg, we deduce that the mixing efficiency is 

Si ~ ^1 (T° F'( V ) 2 dr l + j™ G'(rifd V \ ^ Ve 1 ' 2 . (VI.23) 

Figure 5 also shows the scaling of £\ from the numerical solution of (VI.4). The internal layer 
asymptotic approximation in (VI.23) and the numerical solutions differ by approximately 1% at Pe 
= 1000. The (lack of) scaling in r was also confirmed numerically. Interestingly, for this problem 
stirring at ever smaller scales does not further enhance the mixing efficiency on small scales. This is 
apparently because the decrease in gradient variance due to stirring on small scales is compensated 
by the increase in gradient variance in the internal layer. 

The inverse-gradient variance, 

(IV-^I 2 ) = (\V-\f(y)sm(k s x)+g(y) C o S (k s x))\ 2 ), (VI.24) 
requires a slightly more subtle analysis. Expanding f(y) and g(y) in Fourier series, 

oo oo 

f(v) = ^2fn cos(nk u y), g(y) = ^ g n s\n{nk u y) (VI.25) 
n=0 n=l 
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the inverse gradient variance is 
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(VI.26) 



We know that £_i < Pe, so if we could establish a Zower bound with the same scaling then we 
could conclude that £_i scales classically ~ Pe. Hence we focus on deducing an asymptotic upper 
bound on (jV -1 ^ 2 ). As will be seen, however, the best we can do is assert the classical scaling 
with a logarithmic correction. 

Toward that end, using the fact that Y^=i n ~ 2 = 7I " 2 /6 J 
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(VI.27) 



Then the key is to note that the Fourier coefficients f n are all 0(S/Uk s ) — without any further 
factors of 8 appearing — while the g n are 0{S/Uk s x log 3). Indeed, 
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(VI.29) 



Uk s Jo Uk s 
In the last step above we used the fact that |G(r/)| ~ r]^ 1 as r\ — > 00. 

The unstirred large-scale variance is (|V _1 6>o| 2 ) = S 2 /n 2 k & s , so we deduce the large-scale mixing 
efficiency obeys 

1 



Pe 



(VI.30) 



1 + C" 



log Pc 



where C and C are absolute constants. The numerical solutions of the steady advection-diffusion 
equation shown in Figure 5 confirms this classical scaling; the logarithmic term is not visible for 
this range of Pe. Again, this is consistent with the dynamic data for S-\ in Figure 2. 



VII. SUMMARY AND DISCUSSION 



The suppression of bulk variance of a scalar field is a natural gauge of the efficiency of stirring, 
and this notion allows for the examination of the effect of a flow on various length scales. We have 
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quantified the influence of stirring on weighted bulk variances in terms of nondimensional mixing 
efficiencies £ v that characterize fluctuations on relatively small (p = 1), intermediate (p = 0) 
and large (p = — 1) length scales. We studied these mixing efficiencies for statistically stationary, 
homogeneous and isotropic flow fields stirring scalars sustained by a variety of steady but spatially 
inhomogeneous scalar sources and sinks on periodic domains. We reach a number of conclusions: 

• Very generally, the mixing efficiencies £ v are bounded from below ~ Pe° at high Peclet 
numbers — indicating ineffective stirring — and from above ~ Pe 1 — the classical scaling antic- 
ipated by the simplest eddy diffusion theory. Classical scaling of the efficiencies corresponds 
to the existence of residual supression of variance in the vanishing diffusion (k — > 0) limit. 

• Source, sink and statistically stationary homogeneous and isotropic flow combinations exist 
that realize both the upper and lower bounds' scaling with respect to Pe — and in some cases 
the prefactors as well — showing that the analysis is sharp at the most general level. 

• The small and intermediate scale mixing efficiencies £\ and £q are ultimately limited by 
length scales in the source-sink distributions. That is, even if the efficiencies scale classically 
the inferred mixing lengths ~ nPe/U are limited by length scales defined by the sources and 
sinks rather than by the spectrum of scales in the flow. On the other hand the upper bound 
analysis does not prevent small scales in the flow field from enhancing suppression of the 
variance at large scales, i.e., £-i, to an unlimited degree. 

• Sufficiently "rough" scalar sources and sinks, in particular non-square-integrable distribu- 
tions, necessarily produce sub-classical scaling of £$ and/or £\ at high Peclet numbers. In 
such cases there is no residual suppression of variance in the vanishing diffusion limit for any 
finite mean kinetic energy statistically stationary homogeneous and isotropic stirring. 

• The mixing efficiencies at various length scales need not scale the same as Pe — ► oo. This 
was illustrated explicitly by the example of a monochromatic source distribution stirred by 
the random sine flow. We found by direct numerical simulation and an adiabatic asymptotic 
internal layer analysis that for small scales £\ ~ Pe 1 / 2 , for intermediate scales £q ~ Pe 5 / 6 , 
and the large length scale efficiency £_i exhibits near-classical scaling ~ Pe/logPe. 

The analysis, simulations and modeling presented in this paper have addressed some fundamen- 
tal questions regarding suppression of the long time averaged bulk scalar variance via stirring by 
incompressible and statistically homogeneous and isotropic flow fields, but compelling challenges 
remain for future study. Among the unsolved problems open questions are: 
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Can the anomalous mixing efficiency bounds for measure- valued or fractal sources be achieved 
by any statistically stationary homogeneous and isotropic flow? This is a difficult problem for 
direct numerical simulations of the advection-diffusion equation because of the small spatial 
scales that need to be resolved. An alternate simulation approach is particle tracking, by 
following the evolution of many discrete particles moving with the flow and diffusing and 
approximating the continuous scalar concentration on an appropriately coarse grained level. 

Can the small-flow-scale enhancement of the large-scale mixing efficiency £_i, as suggested 
by the upper bound in (III. 43), be realized? It makes sense that stirring may be capable 
of suppressing large-scale variance, even for negligible molecular diffusion, by transferring 
scalar inhomogeneities from large length scales to small length scales via kinetic stretching 
and folding mechanisms. The upper bound provides a quantitative estimate of this effect 
and it will be interesting to see to what extent it is an accurate estimate. 

What is the mixing efficiency of "real" stationary homogeneous and isotropic turbulence 
in two and three dimensions? The bounds derived here all apply to these flows, and the 
inevitable limitations on the mixing efficiency imposed by the source and sink structure 
apply. Turbulent mixing by a fluid with viscosity v is often characterized by a Reynolds 
number Re = UL/u and a Schmidt number Sc = u/k. While the mixing efficiencies £ p 
are studied here as functions of Pe = Re x Sc, the separate Reynolds and Schmidt number 
dependences are of interest, too. These questions can be investigated theoretically and via 
direct numerical simulations in the periodic domain setting utilized here. 

The uniform-wind optimal solution presented in Section IV is only appropriate for a one- 
dimensional source distribution in a domain with periodic boundary conditions. It is optimal 
by virtue of always having a component blowing along the source gradient, so no time and 
little energy is wasted blowing along lines of constant source. If the source has a more 
complicated structure, then it is not generally possible to achieve such an efficient flow 
while preserving incompressibility. Physical boundaries or a different golbal topology (e.g. a 
spherical domain as relevant to geophysical applications) will likewise invalidate the uniform- 
wind solution. The general question is then: For a given distribution of sources and sinks 
in a given domain, what is the optimal stirring strategy to achieve the greatest reduction 
of bulk variance? This question has obvious relevance to engineering applications. While 
we have answered it in the simplest setting of a single mode source-sink distribution on a 
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periodic domain, it is clearly a complex problem in general. It is not unlikely, however, that 
intuition gained from the study of idealized models may contribute to the development of 
valuable intuition for such systems. 

• Bounded domains with rigid walls are more appropriate for some engineering applications. 
Then the scalar sources and/or sinks can be implemented by boundary conditions rather 
than as body sources and sinks [55]. We can also imagine situations where fluid and/or 
imposed scalar fluxes at boundaries contribute to the scalar variance in the bulk. How this 
will affect the efficiency scalings is an open question. 

• We have seen that very generally S\ > 1, but have only established a similar lower bound on 
So and £_i for the simplest case of monchromatic sources. It is an open question whether 
or not there exist source, sink and flow (even steady flow) combinations where £q and/or 
£_i < 1. It is worthwhile noting in this context that in all cases (6> 2 ) 2 < (iV^l^dV -1 ^! 2 ). 
Any stirring decreases the mean bulk gradient variance, but this suggests that if it does not 
decrease the variance significantly then it would have to increase the inverse gradient vari- 
ance. Hence it may not be surprising to find flows that can amplify large-scale fluctuations. 

• A related problem of interest is the advection-diffusion of a passive scalar 9 with decay rate 
£ sustained by a body source with evolution described by 

f)9 

— +u-V9 = kA9 + s(x)-(9. (VII.l) 

In applications £ may have an interpretation as a chemical reaction rate or a radiative 
relaxation rate relevant in meteorology or for other geophysical flows on the sphere. A linear 
amplitude damping term in the advection-diffusion equation introduces new features such 
as competition between scalar decay and mixing to supress steady-state variances [50]. 

• Finally, all the questions above can be posed for time dependent sources and sinks. 
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APPENDIX A: LOWER BOUND ON VARIANCE EFFICIENCY 

One might expect the variance efficiency £q to have a lower bound of unity, implying that 
stirring always decreases the variance. In order to derive a uniform (in Pe) lower estimate we can 
search for the maximum variance subject to the steady-state production-dissipation constraint: 

(9 2 ) > min {(tf 2 ) | k(|W| 2 ) = («?)}. (A.l) 

The solution of this optimization problem in Fourier space is 

where \x is the Lagrange multiplier enforcing the constraint 

- «?0w = 2 ?^fi- (A - 3) 

where * denotes the complex conjugate. 

In the case of a dichromatic source with wavenumber k\ at amplitude s\ and ki at amplitude 
«2j the constraint requires one to solve a cubic equation for £ = 1/ jink\, 

(1 + a)f + i(l + a/3 + 4/? + 4a)£ 2 + (J3 + a(3 + 1 + a)£ + \{(3 2 + a/3) = (A.4) 

where a = sf/s 2 , and /3 = kf/k 2 . Then the mixing efficiency bound is 

2 > S^I^WI' _ 4(1 + fr) 

and a numerical evaluation of the roots reveals that the minimum value of the efficiency bound is 
less than 1 for V £. Thus the variance production-dissipation balance alone does not rule out the 
existence of flows that could increase the scalar variance. 

APPENDIX B: STEADY UNIFORM WIND ON A 5-FUNCTION SOURCE 

Taking the Fourier transform of the steady advection-diffusion equation (1.3) with a uniform 
velocity field along the x^-axis, 

§(k) = i2 S{k \ Ti (B.l) 
nk 2 + zUkd 
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where is the dth component of the horizontal wavenumber. Now substitute s(k) = const for 
a (5-function source and after approximating the sums by integrals (since we are only interested in 
the asymptotic behavior), we find 

k 2 P +1 dk 



/•27T POO 

d = 2: (\V P 9\ 2 ) = / # / - 
JO JI-kIL k 



d = 3: (\V p e\ 2 ) 



27T/L ^k^ + U 2 ^ COS 2 0' 

2 r: - . f°° k 2 (p +1 Uk 

J 2n/L K 2 k* + U 2 k 2 COS 2 



(B.2a) 
(B.2b) 



The variances in the absence of stirring are found by calculating the integrals (B.2) with U = 0. 
Straightforward evaluation of the integrals (B.2) yields 

d = 2: E x = 1, «S ~v / Pe, £_i ~ VPe, 



d = 3 : £\ = 1, So 



(B.3a) 
(B.3b) 



\/logPe 

These anomalous scalings in Pe suggest that the uniform flow is not the optimal allowed by the 
bound for the 5-function source in both d = 2 and 3. This emphasizes the source-dependent nature 
of the optimal stirrer. 



(C.la) 
(C.lb) 



APPENDIX C: STEADY SHEAR AT AN ANGLE 

When the source is tilted at a 45° angle, i.e. s(x) = V2S sm(k s (x + y)), the functions / and g 
must satisfy 

r d 2 1 

-V2Uk s sin(k u y) g(y) = k -k 2 + f(y) + V2S sin(k s y) , 

r d 2 1 

V2Uk s sm(k u y)f(y) = n -k 2 s + g(y) + V2S cos(k s y) . 
Upon changing variables as done in section VII we obtain 

f(y) + siny = - sm(ry)g(y), (C.2a) 

g(y) + cos y = sm(ry)f(y). (C.2b) 

The internal layer solution is once again obtained by expanding in powers of Ve~ l . To leading 
order the solution to (C.2) in the outer region is 

1 



1 


' 1+^' 

dy 2 


vl 


1 


" 1 

dy 2 


V~e 



/outer ~ 1, pouter ~ ~ tan(ry) ' 

The inner internal layer solution is the same as for the untilted source problem. In Figure 6 we 
compare the internal layer solution against the exact (numerical) solution of the full advection- 
diffusion equation (VI. 4) at Pe = 1000 where the small parameter e w .2. Recovering all the 



0.05 0.1 0.15 0.2 °'° 25 0.05 0.1 0.15 0.2 

y y 

FIG. 6: Comparison of the direct numerical solution (solid) and the internal layer solution (dashed) for 
Pc = 1000. 

scalings and letting 5 = e/r l ^k s , we define the composite approximations 



S S 

/comp(y) = jjj^ [/inner (?/) + /outer ( J/)] = 



k u 5 (<5 ^ 



(C.4a) 



9c °^ {y) = ik Ky x ^ mncr(y) x ° outcAy) = mr. h iM^) G S) • (c - 4b) 

Given the internal layer solution we can compute the multi-scale mixing measures. In fact, to 
leading order the dependence of the efficiencies on the Peclet number is the same. The difference is 
in the prefactors. Figure 7 compares the scalings derived from the internal layer solution (including 
prefactors) to the discrete data from the numerical solution and the numerical solution using the 
random sine flow. 
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